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Cell division: a source of active stress in cellular 
monolayers^ 

Amin Doostmohammadi,^ Sumesh R Thampi,^ Thuan B. Saw,^ Chwee T. Benoit 

Ladoux,^’^ and Julia M. Yeomans*^ 

We introduce the notion of cell division-induced activity and show that the cell division generates 
extensile forces and drives dynamical patterns in cell assemblies. Extending the hydrodynamic 
models of lyotropic active nematics we describe turbulent-like velocity fields that are generated 
by the cell division in a confluent monolayer of cells. We show that the experimentally measured 
flow field of dividing Madin-Darby Canine Kidney (MDCK) cells is reproduced by our modeling 
approach. Division-induced activity acts together with intrinsic activity of the cells in extensile and 
contractile cell assemblies to change the flow and director patterns and the density of topological 
defects. Finally we model the evolution of the boundary of a cellular colony and compare the 
fingering instabilities induced by cell division to experimental observations on the expansion of 
MDCK cell cultures. 

1 Introduction 

The collective migration of cells plays a pivotal function in vi¬ 
tal physiological processes such as embryonic morphogenesisl^ 
tissue developmentl^, wound healingEEj, and tumor growthEEI. 

The collective motion of cells is often explained by considering a 
delicate interplay between biochemical signaling, metabolic pro¬ 
cesses, and mechanical forcesEI. However, the individual function 
played by each of these mechanisms and their relative importance 
remain obscure. 

The sources of mechanical stimuli often originate within the 
cell culture. Examples are cell motility due to molecular mo¬ 
tors deforming the cytoskeleton, and flow fields established by 
cell division. In the past decade, there has been growing evi¬ 
dence of the important impact of mechanical processes such as 
motility-induced forces in the growth and development of cellular 
culturesI^Jtm. The collective migration of cancer cells due to the 
generation of compressive stress during tumor growthEI, the co¬ 
ordination of cell motion by the emergence of a mechanical wave 
during epithelial monolayer expansion!!^, and the emergence of 
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fingering instabilities at the border of a growing tissuel^, are strik¬ 
ing examples. Nevertheless, it is still unclear how the collective 
motion is induced by mechanical stresses generated at the scales 
of individual cells. 

Recently, important correlations have been found between me¬ 
chanical forces within cell cultures and cell proliferation !^^ 1 14 | 
has been shown that the mechanical forces generated by cell di¬ 
vision and apoptosis can result in fluidization of tissue^^ and 
lead to the formation of long-range vortical structures in living 
tissueESl The dependence of cell division rates on internal me¬ 
chanical stress has been suggested as a possible mechanism for 
growth controlEZIand cellular flow generated by proliferation gra¬ 
dients has been shown to be drastically suppressed by mechan¬ 
ical stresses®. Continuum theories have been used to model 
anisotropic tissue growth with active stresses generated by cell 

divisionESESl. 

Despite the widespread implications, the contribu¬ 
tion of cell division to coordinated cell motion remains largely 
unexplored. 

In this paper, we investigate the generation of active stresses 
by cell division. Using generalised equations of nematic liquid 
crystals with the concentration of nematics allowed to vary, we 
reproduce the experimental flow field of dividing Madin-Darby 
Canine Kidney (MDCK) cells and show that a cell division leads 
to a dipole-like flow field at the division site. Even in the absence 
of intrinsic activity of the cells, such division-induced activity can 
generate turbulent-like flows. We study the role of cell division 
on the cellular flow and compare the effects of division-induced 
activity in assemblies of extensile and contractile cells. The role of 
cell division in determining the growth of an interface bounding 
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a layer of cells is also investigated. We find fingering instabilities 
at the cellular front, which are compared to experiments on a 
growing colony of MDCK cells. 

2 Methodology 

2.1 Equations of motion 

There are many models for collective cell migration in the litera- 
t^reHHUl. The cultures of cells are often regarded as a continuous 
medium due to strong contacts between the cellsEEI. Here, in 
order to describe the dynamics of the collective motion of cells, 
we employ a continuum approach based on the hydrodynamic 
equations of active nematicsESEzI 


{dt-\- Ukdk)Qij 3ij — EHij , 

(1) 

3tp + 3i{pUi) = 0, 

(2) 

p{dt-\-Ukdk)ui = djUij, 

(3) 


The molecular field 
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ensures that the system relaxes diffusively to the minimum of a 
free energy, through equation 0 and the diffusion constant F 
controls the time scale over which the relaxation occurs. We use 
the Landau-de Gennes free energy 


^ = \A(q^c-\QijQjif + ^{dkQijf (8) 


where A is a material parameter characterising the coupling be¬ 
tween the nematic order, Q, and the concentration of cells, c and 
K is an elastic constant measuring the energy associated with de¬ 
viations from nematic ordering. In two dimensions there are, in 
general, two elastic constants associated with bend and splay de¬ 
formations, and we take these to be equal. The free energy 0 
corresponds to a molecular field 


where Qij = 2q{ninj — 3ij /2) is the two-dimensional nematic order 
tensor, with director n and magnitude q, taken to characterise the 
coarse-grained orientation of cells, which is a consequence of the 
underlying cytoskeletal network. Based on the comparison of the 
flow field generated in our simulations with that of the experi¬ 
ments, the length scale of a cell corresponds to approximately six 
grid points in the continuum formulation. 

Equation 0 is a convection-diffusion equation for the dynam¬ 
ics of the order parameter field Q. The total density p obeys 
the continuity equation 0- The evolution of the velocity field u 
is described by equation 0, which is the generalisation of the 
Navier-Stokes equations to describe the dynamics of liquid crys¬ 
tals. 

Consider first the evolution of the order parameter field, equa¬ 
tion 0. As a result of the anisotropic shape of the nematogens, 
the director field responds to the gradients in the flow field. This 
is accounted for by the generalised advection term 

Sij ={XEik + COik){Qkj + + {Qik + ^ik/'^){^Ekj — CO^j) 

— 2X{Qij-\-5ij/3){QkidkUi), (4) 

where 

Eij = {diUj T djUf)/2^ COij = {djUi dfUj^/2 (5) 

are the symmetric and antisymmetric parts of the velocity gradi¬ 
ent tensor, respectively called the strain rate and vorticity tensors. 

In equation Q, A is the alignment parameter which deter¬ 
mines whether the collective response of the nematogens to a 
velocity gradient is dominated by the strain or vorticity. The di¬ 
rector aligns at a given angle to a shear flow, called the Leslie 
angle, if 


If this equality is not satisfied the director field rotates (tumbles) 
in a simple shear flow. Generally, A depends on the shape of 
the particles, with A > 0 and A < 0 corresponding to rod-like and 
plate-like particles respectively. 


Hij = -AQij (q^c - (QijQji) /2) + Kid^Qij). (9) 

Turning now to the Navier-Stokes equation the stress term, 
Hij, includes contributions from the viscous stress 

U^i^cous ^2riEij, (10) 

with the viscosity r/, and elastic stresses 

^elastic ^ -P5ij + Qi^Hkj-HikQkj-KdiQudjQkl 

3' 3 ' 

+ A [2{Qij + -^){QkiHik)-Hik{Qkj + 

-{QikA^)iikjl ( 11 ) 

where P is the pressure. The elastic stress is a second consequence 

of the anisotropic nature of the nematogens. It occurs because, as 
the directors turn they induce stresses and hence a contribution 
to the velocity field, often called the ‘back-flow’. In addition, the 
stress term includes the active stress 

( 12 ) 

where the coefficient ^ measures the strength of the activity. In 
general, the active stress includes contributions from stresses in¬ 
duced by molecular motors and actin polymerisation dynamics. 

< 0 describes a contractile system where the flow field of active 
nematogens is along their short axis and C > 0 describes an ex¬ 
tensile system where the flow is along their long axis. The form 
of the active stress tensor is motivated in^^. Here, in order to 
distinguish the division-induced activity from active stresses in 
equation ( p^ , we refer to the latter as ‘intrinsic’ activity of the 
cells. 

Many properties of active nematics can be interpreted by noting 
that the active term appears in the stress, under a derivative, and 
therefore any spatial gradient in the direction or orientation of the 
nematic order induces a flow. An immediate consequence is that 
the active nematic state is hydrodynamically unstableE^^ and 








Fig. 1 Generation of a dipole-like flow field due to a cell division event modelled as a local increase in concentration (left) and experimental 
measurements of the flow field around a dividing cell (right). Red arrows denote velocity fields and the colourmap in the simulation (left) shows 
concentration. The position of the dividing cell in the experiment is marked by the white ellipse. 


the global nematic ordering is replaced by active turbulence, a 
state characterised by swirls and jets in the velocity field and topo¬ 
logical defects in the director fieldl^. Equations. ([^-([^ have 
proven successful in modelling the behaviour of dense suspen¬ 
sions of microswimmers and active suspensions of microtubules 
driven by molecular motorsE^^. 

To introduce cell division, the equation of motion for the con¬ 
centration of cells is taken as 

dfC-\-di{uiC^ = T Ofc, (13) 

where k is the diffusivity of the cells and a represents the prolif¬ 
eration rate due to the growth of cells. We will show that any in¬ 
crease in concentration results in the generation of active stresses, 
which drive the flow of cells. Note that the form of the free energy 
in equation ^ corresponds to a constant concentration through¬ 
out the system in equilibrium. In section 5 of the paper we will 
generalise the free energy and the concentration equation of mo¬ 
tion to allow us to model the mechanical behaviour of a dividing 
cellular colony with a free surface. 

The equations of motion 0 - ([ 3 ]) and ( [13| ) are solved using a 
hybrid lattice Boltzmann methodl^SI. Simulations were performed 
in a two-dimensional square domain of size 200 x 200 and dis¬ 
crete space and time steps were chosen as unity. Unless other¬ 
wise stated, the parameters used in the simulations are p = 40, 
r = 0.1, K* = 0.1, = 0.05, a = 0.0001, A = 0. 3, q=l, C = 0.0025 

and T] = 2/3 in Lattice-Boltzmann unitsE^I^. We prescribe zero 
velocity field, a uniform concentration and random nematogen 
orientations as initial conditions and implement periodic condi¬ 
tions on domain boundaries. 


2.2 Experiments 

Madin-Darby canine kidney (MDCK WT) strain II cells were cul¬ 
tured in low glucose DMEM medium (Invitrogen), with 100 jigjml 
penicillin, 100 fJiglml streptomycin (Invitrogen) and 10% FBS (In¬ 
vitrogen). Experiments for the cell division flow field measure¬ 
ment (DEE) were performed with cells incubated in Leibovitz’s 
L — 15 (Invitrogen), with 100 pg/m/ penicillin, 100 pg/m/ strepto¬ 
mycin (Invitrogen) and 10% FBS (Invitrogen). Under these con¬ 
ditions, the typical aspect ratio of cells (if fitted with an ellipse) 
averaged at ^ 1.8, and could reach up to 3 — 4. 10 pM blebbis- 
tatin (Cayman Chemical Company) drug was added and left in 
the medium throughout the experiment to slow down dynamics 
of tissue and reduce noise in measurements. Measurements were 
done without drug as well, and similar cell-induced vortices were 
also measured, but happening earlier after cell division event. 
The fingering experiment (EE) used a culture medium without 
the drug. 

For DEE, MDCK cells were confined on a 500 pm diameter 
square pattern by a microcontact printing techniquelSI. Confine¬ 
ment is essential to get a confluent tissue without too much spa¬ 
tial variations in cell density, and also to establish nematic order¬ 
ing of cells, since tissue will be polarized if left to expand into 
free area. Fibronectin (FN - 25 pg/m/ Atto dye conjugated FN 
and 100 pg/m/ pure FN, Sigma and Roche) was incubated on 
polydimethylsiloxane PDMS (Sylgard 184; Dow Corning) stamp 
for 1 hr at room temperature, dried and stamped on a PDMS 
spin-coated dish. Pluronics El27 (Sigma) 1% was incubated to 
passivate the non-stamped area, before cells were seeded and left 
overnight to grow to confluence on the stamped area and imag¬ 
ing the next day. For EE, cells were seeded overnight on a FN 
(pure FN, 20 pg/m/, 1 hr incubation) incubated glass bottom dish 







(Ibidi). A PDMS block with long, straight boundaries was placed 
on the surface before cell seeding to prevent cells from penetrat¬ 
ing the desired area. After cells outside the block grew to conflu¬ 
ence, the block was lifted and imaging started when cells started 
to migrate into the void, forming fingers. 

Time-lapse Imaging was done respectively for 2 days (DFF) 
and ^ 10 hrs (FE) (Biostation (Nikon)). Imaging resolution was 
0.647 [imjpixel with a lOX Phase objective. The image interval 
was 10 min. 

The velocity field of a confluent tissue in the confined region 
was measured®Sl. individual cell division events were pin¬ 
pointed by eye. The cell division starting point was chosen to 
be the last frame where the mother cell was still present. Each 
division starting point was chosen to be a new reference point 
(0,0) and the velocity around this point was determined and ro¬ 
tated with respect to the closest nematic director, such that the 
nematic director aligned in the y-axis in this new reference frame. 
The transformed velocity field of division event, /, / time-frames 
after the division starting point is: v/y (the origin of the local 
reference frame at time-frame, /, was always chosen to be the 
same position as the division starting point i.e. at / = 0). The 
average transformed velocity field was determined for different 
/ G 0,1,2,..., as yf(r) = Ef v/j/A. The drift corrected, aver¬ 
age transformed velocity field is y = V/(r) f(r)lNr, 

where Ny is number of local velocity vectors around each new 
origin of area 120 x 120 pn? (one cell length is ~ 15 pm). 

3 Division-induced fiow field and meso- 
scaie turbulence 

The flow field generated by the coordinated motion of cells ex¬ 
hibits turbulent-like flow patterns characterised by a distribution 
of flow vortices, and often referred to as meso-scale or active tur- 
bulencel^3^32l. The emergence of meso-scale turbulence in cellu¬ 
lar matter is induced by the motion of the cells which is driven by 
rearrangements of the cytoskeletal elements by molecular motors. 
This is modeled by the active term, Hij = in the equations 

of motion. In order to separate the effects of the cell division from 
the intrinsic activity of the cells, we first consider the dynamics of 
a concentration of cells without any conventional active forcing 
((^ = 0), but proliferating due to cell division. 

We present results for the flow field around a single dividing 
cell, and then describe the velocity field produced by many divid¬ 
ing cells, at random positions and times, but with a given number 
of divisions, per unit area at any time. Cell division is modelled 
by locally increasing the concentration (by a) within a circular 
area of radius three grid points. The increase in concentration 
is maintained for a short duration of ten lattice Boltzmann time 
steps, much smaller than typical evolution time scales of flows 
in the simulations (^ 10,000). Division events are introduced at 
random positions as a transient, small, local increase in concen¬ 
tration. This locally changes the nematic ordering, because of the 
coupling in the free energy, which in turn drives the flow field. 
We take measurements once the flow has reached a statistically 
steady state. 


3.1 Flow field of a division event 


Recent experimental observations of the flow fields of endothelial 
cells have reported the emergence of velocity vortices around cell 
division sitesCSI. In these experiments the flow field generated by 
the cell division may be affected by a number of mechanical fac¬ 
tors such as the motion of the cells and strong friction with the 
substrate. It has been previously argued that cell division can be 
thought as a local source of energy injectionlSI and stress gener- 
ationl^ and that can be modeled as a force dipolell^. The force 
dipole can be generated as an elongated cell exerts anisotropic 
forces on its neighboring cells. Here we show that the flow field 
is automatically generated in our formulation by a local increase 
in the concentration, which induces a local active stress. This oc¬ 
curs because, as the concentration is increased, the free energy 
drives a corresponding increase in the local nematic order (to¬ 
wards qnem = \/c^) with a time lag of the order tqr^X/T. Though 
the concentration disturbance is isotropic, the nematic symmetry 
of the director field breaks the isotropy. Changes in the local ne¬ 
matic order lead to the variations in the molecular field, Hij, and 
analysis of eqn (11) shows that terms proportional to kQtj domi¬ 
nate in producing an anisotropic elastic stress. 


The flow field generated by a single cell division event is com¬ 
pared with the experimentally measured flow field around a divi¬ 
sion point in MDCK cells (Fig.[^. The flow field reported in Fig.[^ 
is after 200 time steps. In experiments, velocity fields were mea¬ 
sured after 180 min from the cell division and are averaged over 
100 cell division events. The flow field resembles that of a force 
dipole with an octupole correction, which takes into account the 
finite extent of the local increase in the concentration. Moreover, 
as evident from the simulation results in Fig. the division re¬ 
distributes the concentration along the director orientation and 
results in the elongation of the cell parallel to the division axis. 
This is consistent with a recent experiment on MDCK cells, which 
showed that by orientating along the long axis of the cell, division 
redistributes the mass along the long axis and causes the outward 
(inward) movement of neighboring junctions parallel (perpendic¬ 
ular) to the division axisISI. We note that the spatial length scale 
of the flow field generated by a single division event in MDCK 
cells (with reduced motility) spans ^ 100 pm, similar to the ve¬ 
locity correlation lengths measured in MDCK tissue during normal 
proliferation and motilityEni. Taken together, our simple model¬ 
ing approach captures important features of the division event 
and shows that each division acts as a source of activity, generat¬ 
ing hydrodynamic flows and driving the system out of thermody¬ 
namic equilibrium. 


3.2 Meso-scale turbulence of dividing cells 

Figure [^illustrates the flow field of an assembly of dividing cells 
in the absence of any intrinsic activity (i^ = 0). The activity in¬ 
duced by division events drives the system out of equilibrium and 
establishes a spatially and temporally evolving flow field. At a low 
division fraction the flow fields generated by the dividing cells are 
independent, they do not interact each other. (Fig. [^)). At di¬ 
vision fraction of = 0.3, however, the division-induced activity 
results in the formation of persistent regions of high velocities in 
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Fig. 2 Steady state flow fields for a concentration of dividing cells with no intrinsic activity = 0). (a), (b) and (c) correspond to low = 0-05, 
moderate = 0.3, and high = 0.5 division fractions. The flow fields are visualised as streamlines using the Line Integral Convolution (LIC) method, 
the colour maps indicate the magnitude of the velocity and white arrows the local direction of the velocity. We use the same resolution as the flow field 
simulations to construct LIC maps, (c) The rms-velocity of cell assemblies increases with increasing the division fraction (|)^. The sharp increase in 
vrms at (l)d = 0.08 corresponds to the onset of collective behaviour. Panel (d) shows normalised velocity-velocity correlation functions for different 
division fractions. 


straight and curved patterns. The emergence of these jets and 
vortices in the flow field shows that the cell division triggers co¬ 
ordinated patterns of cellular motion (Fig. [^b)). The dynamics 
is reminiscent of ‘active turbulence’, the emergence of collective 
motion in wet active nematics driven by nonzero intrinsic activity 

^ 0)S3 As the division fraction is further increased to = 0.5 
the swirls and jets are still seen, but the length scale of the vor¬ 
tices decreases (Fig.j^c)). 

To quantify these results we measured the root-mean-squared 
(rms) velocity, and the velocity-velocity correlations function as a 
function of the cell division fraction. Results for the rms velocity 
are presented in Fig.j^d). As expected, the rms velocity increases 
as the number of cell divisions increases. Note, however, that 
there is a sharp increase in rms velocity at = 0.08 (red 

vertical line in Fig.j^d)). This is the cell division fraction beyond 
which collective patterns and mesoscale turbulence are observed. 

The transition to mesoscale turbulence is also evident from 
measurements of velocity-velocity correlation function, Cw = 
((v(r,?).v(0,?))/v(0,r)^), where () denotes temporal and spatial 
averaging (Fig.|^e)). Below the critical division fraction the 
correlation length corresponds to that of a single division event. 
Above (j)^, however, the correlation length markedly increases cor¬ 
responding to the emergence of collective behaviour. Further 
increases in result in a gradual reduction of the correlation 
length as the effective activity of the cell assembly increases and 
the jets and swirls shrink in size. 


Although, several mechanical processes such as a bimodal dis¬ 
tribution of cell behavioiESI, compressive pressure® and plitho- 
taxis (cells actively generating heterogeneous intercellular forces 
and migrating in the direction of maximal principal stressESEll), 
have been introduced as potential mechanisms of the collective 
behaviour of cells, the underpinning mechanical source of this 
behaviour is not yet clear. Here we suggest that even cell division 
alone is sufficient to produce stresses, which drive the collective 
migration of cells. 

4 Comparing cell division and intrinsic ac¬ 
tivity 

In our discussions so far, the system operates without extensile or 
contractile stresses due to the intrinsic activity of the cells and the 
sole source of activity is the cell division. In order to characterise 
the combined effects of cell division and the intrinsic activity, we 
next investigate the impact of cell division on the dynamics of an 
assembly of active cells ((^ 0). The underpinning mechanism 

of cell activity is not completely clear, but the prevalent micro¬ 
scopic picture is that the activity is generated through the stress 
that is exerted by myosin motors perturbing the actin cytoskele- 
ton and by the polymerisation of the actin filaments®. Whether 
the stress is extensile or contractile for all type of cells is not yet 
known, but a number of studies show tha t the m yosin contractil¬ 
ity is an effective mechanism of propulsionSSEzl. Moreover, recent 
theoretical predictions suggest that the cell activity can show ex- 
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Fig. 3 The effect of cell division on the dynamics of contractile and extensile cells, (a)-(d), colourmaps show the velocity fields in contractile ((a),(b)) 
and extensile ((c),(d)) systems (visualisation method is the same as in Fig.[^. The colour maps indicate the magnitude of the velocity and white 
arrows the local direction of the velocity, (e)-(h), colourmaps of the director field indicating topological defects (with +1/2 and -1/2 defects denoted by 
red circles and yellow triangles, respectively). (i),(j), effect of cell division on velocity correlation functions for contractile and extensile cells. 













tensile or contractile behaviour depending on the strength of the 
actin—myosin interactionsEH. 

A comparison of the flow fields of cells with and without divi¬ 
sion, shows a notable difference between extensile and contractile 
cells (Fig. 1^. Figure [^a) shows the active turbulent behaviour 
induced by contractile activity of the cells < 0). It is compared 
with Fig. [^b), where cell division events are incorporated. As 
discussed earlier, the cell division events are associated with ex¬ 
tensile stresses. The extensile and contractile contributions act in 
opposition to reduce the net stress in the system. The rms ve¬ 
locity falls by ^ 70% corresponding to a much reduced effective 
activity. In contrast, introducing division to a culture of extensile 
cells, the two contributions to the activity add up resulting in an 
increase in the stress, and the rms velocity inceases by ~ 125% 
(Fig.|5c),(d)). 

The cell division locally generates vortical structures. These 
smaller vortices thus appear inside the larger swirls generated 
by the effective activity, producing an assembly that is charac¬ 
terised by large jets and swirls, which are interleaved by small 
vortices generated by cell division. This is a characteristically dif¬ 
ferent flow pattern than that observed in meso-scale turbulence 
due purely to the intrinsic activityESI. 

We also calculate the correlation functions for velocity fields 
(Fig. I^i), (]■)). The effect of cell division on the velocity cor¬ 
relation length shows opposite trends in contractile and exten¬ 
sile assemblies: the correlation length is increased in the former 
(Fig. [^i)), while it is reduced in the latter (Fig. This fol¬ 

lows the same trend as recent results for extensile active nemat¬ 
ics, driven by the intrinsic activity, in the absence of an ordering 
free energy which found that the correlation length increases (de¬ 
creases) with decreasing (increasing) activityl^. However, here 
the variation of the cell concentration due to division events in¬ 
troduces new dynamics to the system and any analogy must be 
treated with caution. 

Another important consequence of cell division is in changing 
the number density of topological defects, which may be impor¬ 
tant in controlling the structure of cell layers. Recent experiments 
on fibroblasts cells have shown that the nematic order of the cells 
is accompanied by formation of topological defects, which pre¬ 
vent the development of infinite size nematic domainsEU. In the 
simulations, the number density of topological defects increases 
with increase in cell division in a system with no intrinsic activity 
((^ = 0). However, the cell division has a different effect when 
it is associated with systems having intrinsic extensile or con¬ 
tractile activities. While the generation of topological defects is 
enhanced by cell division in extensile systems (Fig. [^g),(h)), it 
is significantly reduced in contractile assemblies (Fig. |^e),(f)). 
This again can be explained by the division reducing (increasing) 
the effective activity of contractile (extensile) active nematics in 
accord with recent studies showing that the number of defects in¬ 
creases with increase in the activity of the systemEH. Taking the 
effects of cell division, a decrease (increase) in the number of de¬ 
fects for contractile (extensile) systems is observed since effective 
activity is reduced (increased). Less topological defects implies 
less stress in the tissue, which might have important physiologi¬ 
cal implications. 
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Fig. 4 (a) Temporal evolution of a free surface of MDCK cells and 
emergence of the fingering instability in experiments (left) and the same 
phenomena observed in our simulations by division-induced activity 
(right). In the simulations, colormaps show the concentration. The time 
step in experiments is 150 min and in the simulation it is 120 in simulation 
units, (b), (c). A close-up of the velocity field in the band for experiment 
and simulation, respectively. 


5 Cell division and the free surface 


Up to now, we have considered division effects on the dynam¬ 
ics of cell assemblies in periodic domains. In many physiological 
applications such as morphogenesis, tissue expansion, and wound 
healing, the mechanical response of a free surface to the cell inva¬ 
sion is of considerable importanceEHZI. Here, using the equations 
of lyotropic active nematics, we extend our results to the case 
where a cell assembly is separated from an otherwise isotropic 
liquid by a free interface. To distinguish the cell culture from the 
isotropic fluid, we define a scalar order parameter 0, which mea¬ 
sures the relative density of each component with 0 = 1 for the 
cells and 0 = 0 for the isotropic fluid and evolves according to the 
Cahn-Hilliard equationEH 




( 14 ) 







































where is the mobility, = 5^/50 is the chemical potential 
and the free energy of the system is 


+ 2 2^{^kQij)^: 


(15) 


where A 0 and are material constants. Equation (14 ) together 


with equations 0-0 are solved here to describe the dynamics 
of a dividing colony of cells with free surfaces. An additional term 
Uij = (^ - /d{dj^)) must be added to the stress 

components in eqn. when the variable (j) is introduced. More 
details of the form of the free energy and the governing equations 
of lyotropic active nematics can be found inE3. We use r 0 = 0 . 1 , 
A(j) = 0.08 and Kfj, = 0.01. We do not explicitely include any terms 
in the free energy that lead to interface anchoringEj, but active 
anchoring may result from hydrodynamic stresses at the inter- 
faceESl. 

In Fig.ga), numerical results for the time evolution of the sur¬ 
face of a cellular layer are compared to the results of experiments 
on the growth of the surface in a colony of dividing MDCK cells. 
Unlike the experiments, we consider cells with no intrinsic ac¬ 
tivity ((J = 0 ) in the simulation to show that a similar behaviour 
follows from considering the division-induced activity alone. Pre¬ 
vious studies have predicted that existence of source terms such 
as material production can drive hydrodynamic instabilities in the 
form of undulations at the interface between a viscous fluid and 
viscoelastic materialES. As evident from Fig.j^a), the expansion 
of the band is accompanied by instabilities that lead to the forma¬ 
tion of fingers at the surface in both experiment and simulation. 
Although previous studies have associated the fingering instabil¬ 
ities to the formation of leader cells at the bordei®, our results 
suggest that the same phenomena can be induced due to the in¬ 
stability of the nematic field to division-induced activity. It is well 
known that the presence of activity can result in the formation 
of bend instabilities in extensile active nematicsE^ESI. Since the 
cell division introduces extensile stresses to the cell culture, it 
can lead to the instability of the nematic field of the cells and 
induce instabilities at the surface. In addition to the emergence 
of fingering instabilities at the surface, long-range velocity fields 
are generated within the growing band even far away from the 
surface (Fig. [^b),(c)). The appearance of long-range velocity 
fields with no preferred orientation towards the free surface has 
been reported in previous studies of tissue growth in response to 
a model wound^^. However, in explaining the experimental ob¬ 
servations, the emergence of long-range velocity fields and their 
correlations with the cell movements were attributed to ‘cryptic’ 
lamelliopodia, which spread underneath other cells during the tis¬ 
sue growthlS^, while here the collective motion is induced by cell 
division. Taken together, the comparison of our simulations with 
experimental observations on the dynamic evolution of the free 
surface of a cell culture shows that similar qualitative behaviour 
such as fingering instabilities at the border and long-range veloc¬ 
ity fields can be induced by the cell division-induced activity. 


6 Conclusions 

To conclude, we propose a modeling framework that describes the 
effect of cell division on the dynamics of cell cultures and demon¬ 
strate that the model reproduces the experimentally measured 
flow field around a dividing cell in a MDCK cell culture. We show 
that an extensile active stress can naturally arise from the nema- 
tohydrodynamic representation of cells due to a local increase in 
cell concentration. We demonstrate that even in the absence of 
active forcing due to intrinsic activity of the cells, the cell divi¬ 
sion alone can lead to a coordinated motion of cells. The results 
suggest that cell division can be considered as one regulator of 
activity in cultures of extensile and contractile cells. Moreover, 
we show that the dynamic evolution of a free interface due to 
the division-induced activity alone (without the intrinsic activity) 
resembles the experimental observations of the expansion of the 
cells resulting in fingering instabilities at the interface, as well as 
previous results on tissue growth and wound healingl^. 

The concept of division-induced activity leads to a broad range 
of questions about the mechanics of growth in cellular assemblies, 
for example the escape of cellular layers into the third dimension. 
The emergence of collective motion due to the active stress gener¬ 
ated by cell division alone can stimulate new mechanisms for the 
control and guidance of cell migration. Our predictions suggest 
experiments aimed at elucidating the mechanical effects of cell 
division on salient features of cell assemblies such as the emer¬ 
gence of glassy behaviour at large concentration of cells and the 
propagation of waves during tissue growth. Future studies will 
focus on a direct comparison of the simulations and experiments 
to further characterize the division effect on flow structure and 
topological defects. For example, a recent study has shown that 
the frictional damping of the momentum on the scale of hydro- 
dynamic screening length, ^ y^q/py, set by the competition of 
viscosity rj and friction coefficient y, can play an important role 
in the dynamics and pattern formation in active nematicsES. The 
substrate friction will be an important contribution to the dynam¬ 
ics of cellular layers that we hope to investigate in the future. 
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